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Abstract. Motivated by emerging applications from imaging processing, the heat flow of a 
generalized p-harmonic map into spheres is studied for the whole spectrum, 1 < p < oo, in a unified 
framework. The existence of global weak solutions is established for the flow using the energy method 
together with a regularization and a penalization technique. In particular, a BV-solution concept 
is introduced and the existence of such a solution is proved for the 1-harmonic map heat flow. The 
main idea used to develop such a theory is to exploit the properties of measures of the forms A ■ Vv 
and A A Vv; which pair a divergence-L 1 , or a divergence-measure, tensor field A, and a BV- vector 
field v. Based on these analytical results, a practical fully discrete finite element method is then 
proposed for approximating weak solutions of the p-harmonic map heat flow, and the convergence of 
the proposed numerical method is also established. 
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1. Introduction. Let fl C R m be a bounded domain with smooth boundary 
<9f2, and S 1 " -1 denote the unit sphere in R n . A map u £ C 1 (ri,S'" _1 ) is called a 
p-harmonic map if it is a critical point of the following p-energy 

EJv) :=- [ \S7v\ p dx, l<p<oo. (1.1) 
P Ju 

It is well-known, [T2J [2TJ [3D] , that the Euler-Lagrange equation of the p-energy is 

-A p u=|Vu| p u, (1.2) 

where 

A p u :=div(|Vu| p - 2 Vu). (1.3) 



Note that A p is often called the p-Laplacian. It is easy to see that equation (1.2) is a 
degenerate elliptic equation for p > 2, and a singular elliptic equation for 1 < p < 2. 
These degeneracy and singular characteristics both disappear when p = 2. 

We call a map u £ W 1,p (£l, S 11 ^ 1 ) a weakly p-harmonic map if u satisfies equation 
(1.2 1 in the distributional sense. Here W ,p (fl, S n ) denotes the Sobolev space 

W^fi.S"- 1 ) := {u e W l ' p (Sl, R n ); u(x) £ S" 1 " 1 for a.e. x £ 0} . 

One well-known method for looking for a weakly p-harmonic map is the homotopy 
method (or the gradient descent method) , which then leads to considering the follow- 
ing gradient flow (or heat flow) for the p-energy functional E p 

u t - A p u= |Vu| p u in Cl T := Q x (0,T), (1.4) 
|u| = 1 in Or, (1.5) 
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complemented with some given boundary and initial conditions. Clearly, equation 
(1.4) is a degenerate parabolic equation for p > 2, and a singular parabolic equation 
for 1 < p < 2. Again, these degeneracy and singular characteristics both disappear 
when p = 2. 

We remark that p-harmonic maps and weakly p-harmonic maps between two 
Riemannian manifolds (M,g) and (N,h), and their heat flows can be defined in the 
same fashion (cf. [5T][35]). In this paper, we shall only consider the case M = f2 and 
N = S'™ -1 , which is sufficient for the applications that we are interested in. 

The p-harmonic map and its heat flow, in particular, the harmonic map and the 
harmonic map heat flow (p = 2), have been extensively studied in the past twenty 
years for 1 < p < oo (cf. [TTJ H2 H3 EEJ HHJ ETJ H2 |30l |32j [31j [34l [Ml IM1 SQl EJ |46l ^ 
for 1 < p < oo; [55J [57J for p = 1). In the case when the target manifold is a sphere, 
the existence of a global weak solution for the harmonic flow was first proved by Chen 
in |llj using a penalization technique. The result and the penalization technique 
were extended to the p-harmonic flow for p > 2 by Chen, Hong and Hungerbiihler 
in [12]. The p-harmonic flow for 1 < p < 2 was solved by Misawa in [41] using a 
time discretization technique (the method of Rothe) proposed in [31], and by Liu 
in |39] using a penalization technique similar to that of |12j . The p-harmonic flow 
(1 < p < oo) from a unit ball in R m into S 1 C R 2 was studied by Courilleau 
and Dcmcngcl in [16 . In the case of general target manifolds the p-harmonic flow 
(1 < p < oo) with small initial data was treated by Fardoun and Regbaoui in [22], and 
the conformal case of the p-harmonic flow was considered by Hungerbiihler in [33 . 
Nonuniqueness of the p-harmonic flow was addressed in |151 132j . Recently, Giga et 
al. [55] proved existence of strong local solutions for 1-harmonic map heat flow using 
nonlinear semigroup theory. Besides the great amount of mathematical interests in 
the p-harmonic map and the p-harmonic flow, research on these problems has been 
strongly motivated by applications of the harmonic map and its heat flow in liquid 
crystals and micromagnetism, we refer to [TJ [5J EH HHJ (Ml 133 US] and the references 
therein for detailed expositions in this direction. 

Another reason, which is the main motivation of this paper, for studying the 
p-harmonic map and its heat flow, in particular, for 1 < p < 2, arises from their 
emerging and intriguing applications to image processing for denoising color images 
(cf. [49j [50]). We recall that a color image is often expressed by the RGB color 
system, in which a vector I(x) = (r(x), g(x), b(x)) for each pixel x = (2:1, Xg) is used 
to represent the intensity of the three primary colors (red, green, and blue). The 
chromaticity and brightness of a color image are deduced from the RGB system by 
decomposing I(a;) into two components 



77(2) := I I(.t) I, (brightness) 

I(aO , _ . . . 

g := , chromaticity 

\ l K x ) I 

where |I(x)| stands for the Euclidean norm of l(x). By definition, the chromaticity 
must lie on the unit sphere S 2 . One of the main benefits of the chromaticity and 
brightness decomposition is that it allows one to denoise rj and g separately using 
different methods. For example, one may denoise r\ by the well-known total variation 
(TV) model of Rudin-Osher-Fatemi [15] (also see [53]), but denoise g by another 
model. One such model is to define the recovered chromaticity u as a (generalized) 
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p-harmonic map [49j [50 



u= argmin J p ,a( v ) for p > 1, (1.6) 



where 



J PiA (v) := E p (v) + £ / |v - g| 2 for A > 0. (1.7) 

1 Jn 

In particular, the cases 1 < p < 2 are the most important and interesting since the 
recovered images keep geometric information such as edges and corners of the noisy 



color images. We shall call (1.6 1 the p-harmonic model for color image denoising. 



We also remark that the second term on the right-hand side of (1.7 1 is often called 
a fidelity term. As in the TV model [HI [22], the parameter A controls the trade-off 
between goodness of fit-to-the-data and variability in u. 



Again, to find a solution for the p-harmonic map model (1.6 1, we consider the 



gradient flow (heat flow) for the energy functional J Pi \, which is given by 

u t - A p u + A(u - g) = fi P! \u mCl T , (1.8) 

|u| = 1 inft T , (1.9) 

B p n = on dfl T := dfl x (0,T), (1.10) 

u = u onftx{t = 0}, (1.11) 

where 

£ p :=|Vur 2 Vu, Mj) , A :=|Vu| p + A(l-u-g). (1.12) 

The goals of this paper are twofold. Firstly, we shall present a general theory of 
weak solutions for the parabolic system ( | 1 .8[ )- ( p~TTT] ) for the whole spectrum 1 < p < 
oo. To the best of our knowledge, there is no theory known in the literature for the 
1-harmonic map and its heat flow, which on the other hand is the most important 
(and most difficult) case for the color image denoising application. So our theory and 
result fill this void. Furthermore, our theory handles the p-harmonic map heat flow 



( |1.8| -( 1.11 1 for all 1 < p < oo in a unified fashion, rather than treating the system 
separately for different values of p and using different methods (cf. [TTJ [12j [13j [16j UHl 
EH [22J [301 EU E4l SU [461 148]). Secondly, based on the above theoretical work, we also 
develop and analyze a practical fully discrete finite element method for approximating 
the solutions of the p-harmonic map heat flow ( |1 . 8 [ ) - ( p~TTT| ) . 

We now highlight the main ideas and key steps of our approach. Notice that 
there are two nonlinear terms in the p-harmonic flow: the p-Laplace term and the 
right-hand side due to the nonconvex constraint |u| = 1, so the main difficulties are 
how to handle these two terms and how to pass to the limit in these two terms when a 
compactness argument is employed. To handle the degeneracy of the p-Laplace term, 
we approximate the p-energy E p (y) by the following regularized energy 

^(v):=M^ / |Vv| 2 dx+- / |Vv|^ (1.13) 
L Jn P Jn 

= /j^|Vv| 2 + i[|Vv| 2 + a p (e) 2 ] f }d X , 
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where e > and 
Op(e) := 



if 2 < p < oo, 
if 1 < p < 2, 



6 p (e) := e a for 1 < p < oo 



(1.14) 



for some a > 0. Here and in the rest of this paper we adopt the shorthand notation 

|Vv| £ := ,/|Vv|2 + ap (e)2. (1.15) 



To handle the nonconvex constraint |u| = 1, we approximate it by the well- 
known Ginzburg-Landau penalization [7], that is, we abandon the exact constraint, 
but enforce it approximately by adding a penalization term to the regularized p-energy 
Ep. To this end, we replace the energy E e p by 



E E /(v) :=££(v) + L a (v) fore,5>0, 



(1.16) 



where 



L*(v) 



F(v) dx, 



F(v) := 7 (|v| 2 -i; 



V5>0, veR™. (1.17) 



So the idea is, as (5 gets smaller and smaller, the energy functional E^' S becomes more 
and more favorable for maps u which take values close to the unit sphere 

Consequently, the regularized model for the p-harmonic model ( 1.6 1 (with general 
m and n) reads 



where 



u e,<5 = argmin J £, J(v) for p > 1, 



(1.18) 



(1.19) 



In addition, the gradient flow for the regularized energy functional J^'f is given by 



U, 



p 



e,6 i 2_ l)u e,5 



A(u e ' 5 - g) = 



B £ /n = 



in fi T , (1.20) 

on dfl T , (1.21) 

u onfix{t = 0}, (1.22) 

which is an approximation to the original flow ( |1.8| -( |L1~1~| ), where 



u 



e,S 



Bf := 6 p (e)Vu' 



|Vu 



e,5ip-2 



Vu E 



A^u^" := divB p 



■e,S 



(1.23) 



After having introduced the regularized flow ( 1.20 )-( 1.22 ), the next step is to 
analyze this regularized flow, in particular, to derive uniform (in e and 5) a priori 
estimates. Finally, we pass to the limit in ( 1.20 1-( 1.22 1, first letting 5^0 and then 
setting e — > 0. As pointed out earlier, the main difficulty here is passing to the limit 
in two nonlinear terms on the left-hand side of (1.201. For 1 < p < oo, this will be 
done using a compactness technique and exploiting the symmetries of the unit sphere 
S"- 1 , as done in [TU HH [T31 EH E21 HH US] . However, for p = 1, since L 1 ^) is not a 
reflexive Banach space, instead of working in the Sobolev space W 1 ' 1 ^), we are forced 
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to work in BV(fl), the space of functions of bounded variation, since solutions of the 
original 1-harmonic map heat flow |L8|-( |1.11| | only belong to L°°((0,T); [BV(n)] n ) 
in general. This lack of regularity makes the analysis for p = 1 much more delicate 
than that for 1 < p < oo. 



We note also that the regularized flow ( 1.20 1-( 1.22 1 not only plays an important 
role for proving existence of weak solutions for the flow ( |1 .8[ )-( [TTTT| ) , but also provides 
a practical and convenient formulation for approximating the solutions. The second 
goal of this paper is to develop a practical fully discrete finite clement method for ap- 



proximating solutions of the regularized flow (1.20 1-( 1.22 1; and hence, approximating 
solutions of the original p-harmonic flow ( 1.8 )-( 1.11) via the regularized flow. It is 
well-known that explicitly enforcing the nonconvex constraint |u| = 1 at the discrete 



level is hard to achieve. The penalization used in ( 1.20 1 allows one to get around this 



numerical difficulty at the expense of introducing an additional scale S. As expected, 
the numerical difficulty now is to control the dependence of the regularized solutions 
on 8 and to establish scaling laws which relate the numerical scales (spatial and tem- 
poral mesh sizes) to the penalization scale 8 for both stability and accuracy concerns. 
We refer to [B] and the reference therein for more discussions in this direction and 
discussions on a related problem arising from liquid crystal applications. Borrow- 
ing a terminology from the phase transition of materials science, the regularized flow 
( |1.20 1-( 1.22[ ) may be regarded as a "diffuse interface" model for the original "sharp 
interface" model ( |1.8[ )-( fLlT] ), and the "diffuse interface" is represented by the region 

{|u«. ff | <!-*}; 

The remaining part of the paper is organized as follows. In Section [2] we collect 
some known results and facts, which will be used in the later sections. In Section [3] 
we present a complete analysis for the regularized flow ( 1.20 )-( 1.22 ), which includes 
proving its well-posedness, an energy law, a maximum principle, and uniform (in e 
and 8) a priori estimates. As expected, these uniform a priori estimates serve as the 
basis for carrying out the energy method and the compactness arguments in the later 
sections. In Section E] we pass to the limit in ( 1.20 1-( 1.22 ) as 8 — > for each fixed 
e > 0. As in [T3l 1X21 1321 |4T1 . the main idea is to exploit the monotonicity of the 

m 



[Hi nana e] 

operator Af, and the symmetries of the unit sphere S n ~ 



Sections 



are devoted to 



passing to the limit as e — > in the e-dependent limiting system obtained in Section 
g) For 1 < p < oo, this will be done by following the idea of Section [4j However, for 
p = 1 the analysis becomes much more delicate because the nonreflexivity of VF 1:1 (f2) 
forces us to work in the BV(£l) space. The main idea used to develop a BV solution 
concept is to exploit the properties of measures of the forms A- Vv and .AAVv; which 
pair a divergence--^ 1 , or a divergence- measure, tensor field A, and a -B^-vector field v. 
Finally, based on the theoretical results of Sections [3][6] in Section [7] we propose and 
analyze a practical fully discrete finite element method for approximating solutions of 
the p-harmonic flow ( |1.8| -( [fTT] ) via the regularized flow (1.20 1-( 1.22| . It is proved that 
the proposed numerical scheme satisfies a discrete energy inequality, which mimics the 
differential energy inequality, and this leads to uniform (in e and 8) a priori estimates 
and the convergence of the numerical approximations to the solutions of the flow 
(1.8l-(l.lll as the spatial and temporal mesh sizes, and the parameters e and 8 



all 



tend to zero. 



2. Preliminaries. The standard notation for spaces is adopted in this paper 
(cf. [H[3]). For example, W k,p (£l), k > integer and 1 < p < oo, denotes the Sobolev 
spaces over the domain f2 and || • \\w k <p denotes its norm. W°' p (£l) = L p (£l) and 
W k < 2 {&) = H k (n) are also used. L«((0,T); W k > p (Sl,R n )) denotes the space of vector- 
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valued functions (or maps), whose W ,p (Q)-norm is L 9 -integrable as a function of t 

T 1 

over the interval (0,T), and || • Hl^i^'P) := (J Q \\ • W^k.pdt) 5 , for q S [l,oo), denotes 
its norm; and with the standard modification for q = oo. We use (•, •) to denote the 
standard inner product in L 2 (fl), and (■, ■) to denote a generic dual product between 
elements of a Banach space X and its dual space X' . 

In addition, A4(f2) (resp. [A1(f2)] n ) denotes the space of real-valued (resp. Re- 
valued) finite Radon measures on f2. Recall that A4(£l) is the dual space of Co(f2) (cf. 
[3]). For a positive non-Lebesgue measure /i £ A4(f2), L p (£l,n) is used to denote the 
space of L p -integrable functions with respect to the measure fi, and for / £ L 1 (n,/i), 
the measure ffi is defined by 

fn{A) := / f dfj, for any Borel set A C fl. 
J A 

For any /i £ \i = ji a + fi s denotes its Radon-Nikodym decomposition, where 

\x a and [X s respectively denotes the absolute continuous part and the singular part of 
\i with respect to the Lebesgue measure C n . 

In addition, BV(£l) is used to denote the space of functions of bounded variation. 
Recall that a function u £ L (Q) is called a function of bounded variation if all of 
its first order partial derivatives (in the distributional sense) are measures with finite 
total variations in £1. Hence, the gradient of such a function u, denoted by Du, is a 
bounded vector- valued measure, with the finite total variation 

\Du\ ee \Du\(Q) :=sup{^ -udivvdx; v £ [C^(Q)] n , ||v|| L = < 1 } . (2.1) 

BV(fl) is known to be a Banach space endowed with the norm 

\\u\\ BV :=\\u\\ L x + \Du\. (2.2) 

For any u £ BV(il), (Du) a and |-Du| a are used to denote respectively the absolute 
continuous part of Du and \Du\ with respect to the Lebesgue measure C n . We refer 
to [3] for detailed discussions about the space BV(£l) and properties of BV functions. 

For any vector a £ R™ and any matrices A, B £ R nxm with jth column vectors 
A^ and respectively, we define the following wedge products 



^Aa:=^ (1) Aa,--.^ (m) Aa], a A A := [a A • • • ,a A A [m) ] (2.3) 
A A B := A {1) A + • • • + A {m] A B (m) . (2.4) 



We point out that (2.4 1 defines A A B to be a vector instead of a matrix, which seems 
not to be natural. However, we shall see in Section [6] that it turns out that this is a 
convenient and useful notation. 

We conclude this section by citing some known results which will be used in the 
later sections. The first result is known as " the decisive monotonicity trick", its proof 
can be found in [52"] . 

Lemma 2.1. Let X be a reflexive Banach space and X' denote the dual space of 
X . Suppose that an operator T : X — > X' satisfies 

(i) T is monotone on X , that is, (J-u — J-v, u — v) > for all u, v £ X ; 

(ii) T is hemicontinuous, that is, the function t i— > (J-(u + tv),w) is continuous 
on [0, 1] for all u,v,w £ X . 
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In addition, suppose that Uk — > u and Tu k — > / weakly in X and X' , respectively, as 
n — » oo, and 



T/ie 



lim {Tu kl u k ) < (f,u). 



Tu=f. 



The second lemma is a compactness result, it can be proved following the proof 
of Theorem 2.1 of [T5] (also see Theorem 3 of Chapter 4 of [2D]) using the fact that 
the operator A £ is uniformly elliptic. 

Lemma 2.2. For 1 < p < oo, let p* = max{2,p}. For a fixed e > 0, Zei {w e,,5 },5>o 
6e bounded in L°°((0, T); W 1 ^* (0,R™)) and {f e '*}5 >0 &e bounded in L x ((0, T); i x (f2, 
R")), &o£/i uniformly in S. Moreover, suppose that { d ™ t }s>o is bounded in L 2 ((0,T); 
L 2 (f2,R")) uniformly in 5 and w e ' S satisfies the following equation 

dw E > s 



dt 



e,<5 



5>0 



Then {w E ' 5 }«j >0 is precompact in L 9 ((0, T); 14 7l,9 (fi, R™)) 
The last lemma is a variation of Lemma 2.2 its proof can be found in |38j (also 



in the distributional sense 
for all 1 < q < p* . 



see Lemma 9 of [IPJ). 

Lemma 2.3. For 1< p < oo, let {w e } £>0 be bounded in L°°((0, T); LF 1,p (r2, R n )) 
and {f e } &e bounded in -/^((O, T); i 1 (fi, R™)), &o£/i uniformly in e. Moreover, suppose 
that {^-} is bounded uniformly in e in L 2 ((0, T); L 2 (f2, R™)) and w e satisfies the 
following equation 



in fin 



e > 



T/ien {w £ } £>0 precompact in L 9 ((0,T); W 1>9 (ft,R n )) 



in £/ie distributional sense, 
for all 1 < q < p. 

Remark 2.1. In i/ie Lemma \2.2\ e > is a fixed parameter and the differential 
operator A £ does not depend on the variable index S. It can be shown that (cf. Theo- 
rem 2.1 of f!2f) that the lemma still holds for e = when p > 2. On the other hand, 



e is a variable index in Lemma 2.3 and the operator A £ also depends on the variable 
index e. 



3. Well-posedness of the regularized flow ( 1.20 1-( 1.22 1. In this section we 



shall analyze the regularized flow (1.20 1-( 1.22 ) for each fixed pair of positive numbers 



(e, S). We establish an energy law, a maximum principle, uniform (in both e and S) a 
priori estimates, existence and uniqueness of weak and classical solutions. We begin 



with a couple of definitions of solutions to ( 1.20 )-( 1.22 ) 



Definition 3.1. For 1 < p < oo, a map u 



e,S 



tin 



R" 



is called a global weak 
solution to fll.20| -( |1.22| if 

(i) u £ > 5 G J L oo ((0,T);M /1 'P*(fi,R™))n J ff 1 ((0 I 71 ); i2 (^ I R-™)); forp* = max{2,p}, 

(ii) |u e,4 | < 1 a.e. in tt T , 

,e,6 



if 



(iii) u £ ' satisfies ( 1.20 1 -( 1 .22 1 in the distributional sense. 
Definition 3.2. A weak solution u e > s to ( 1.20 )-( 1.22 1 is called a strong solution 
u £ ' s G £ p ((0,T); W 2 >p* (n,R n )) n H 1 ((0,T);LP'(n,H n )). It is called a regular 



solution if in addition u £ ^ s G H X ((Q,T); W 1 * (ft, R")). 
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3.1. Energy law and a priori estimates. Since ( 1.20 )-( 1.22 1 is the gradient 



flow for the functional J^\, its regular solutions must satisfy a dissipative energy law. 
Indeed, we have the following lemma. 

Lemma 3.3. Let Uo and g be sufficiently smooth, and suppose that u e S is a 
regular solution to ( 1.20 1-( 1.22b, then u £ ' 5 satisfies the following energy law 



j;i(u e ' 5 (*)) + I II <'°(t) \\h dt = J$(uo) for a.e. s E [0,T]. (3.1) 



.s,S, 



7 e,8 



Proof. Testing equation (|1.20| with u^' 5 we get 



\<'\t)\\h 



d 
dt 



/{^M|Vu^(t)| 2 + Vu^(t)|? 
Jsi 1 z P 

+ 1 -F(u^(t)) + ^(t)- g \ 2 }dx = 0. 



The desired identity (3.1 ) then follows from integrating the above equation in t over 
the interval [0, s] and using the definition of Jt\- □ 

The above energy law immediately implies the following uniform (in e and S) a 
priori estimates. 

Corollary 3.4. Suppose that Uq and g satisfy 



</p iA (u ) < c 



(3.2) 



for some positive constant Co independent of e and 5, then there exists another positive 
constant C :— C(p,X,Cq) which is also independent of e and S such that 

|| u e ' a || L » (w i,p ) + || u 6 - 4 11^1^2) <C / O rl<p<co, (3.3) 



,e,<5|2 



- 1 llz,~(L 2 ) < C forl<p<oo, 



\Vu s - s \P- 2 Vu s - s \\ Loo(Lpl) <C forp' 



P -r 



1 < p < 00, 



b p (e) || Vu e ' d <C for 1 < p < oo. 



(3.4) 
(3.5) 

(3.6) 



Next, using a test function technique of [12], we show that the modulus \u £ ' S \ of 
every weak solution u e ' S to ( 1.20 1-( 1.22 1 satisfies a maximum principle. 

Lemma 3.5. Suppose that |g| < 1 and \uq\ < 1 a.e. in Q. Then weak solutions 
to (II. 20} 41. 221 satisfy \u 6 ' S \ < 1 a.e. in fi T . 



Proof. Define the function 





z - 1 



for < z < 1, 
for z > 1. 



(3.7) 



It is easy to check that x is a nonnegative monotone increasing function on the interval 
[0,oo). 



Now testing ( |1.20| with v := u s ' S x(\u s ' S \) we get 
1 d 



6\2 



, dx+ f {6 p (e)|Vu^| 2 + |Vu^|r 2 |Vu- , 

l § ( |u^| 2 - 1 )|u^| 2 + A(|u^| 2 - • g) } x(|u^|) dx 
[6p(e) + |Vu E ' g |g- 2 ] |V|u £ - 51212 

{|u«.«|>l} 



u 



e,c5|3 



= 0. 
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Since u e,<5 -g < |u e,(5 |-|g| < |u e,<5 |, the second integral is nonnegative. The assertion 
then follows from integrating the above inequality and using the assumption | u 1 < 1 
a.e. in f2. □ 

3.2. Existence of global weak and classical solutions. We now state and 



prove the existence of global weak and classical solutions to the regularized flow ( 1.20 1- 
(1.221 for each fixed pair of positive numbers (e,S). Since — is uniformly elliptic 



for all 1 < p < oo, the existence of classical solutions follows immediately from the 
classical theory of parabolic partial differential equations (cf . [35 ) . 

Theorem 3.6. Let C R m be a bounded domain with a smooth boundary. 
Suppose that Uq and g are sufficiently smooth functions (say, Uo,g G [C 3 (il)] n ) and 
satisfy (3.2 1. Then for each fixed pair of positive numbers (e, <5), the regularized flow 
( 1.20 1-( 1.22 1 possesses a unique global classical solution u £ ' S . Moreover, u £ ' S satisfies 



the following energy law 



.e,<5 



for a.e. s S [0, Tj. 



(3.8) 



Proof. The existence and uniqueness follow from an application of the standard 
results for parabolic systems, see Chapter 5 of [35] ■ (3.8) follows from Lemma 3.3 □ 



For less regular functions Uo and g, we have the following weaker result. 

Theorem 3.7. Let C R m be a bounded domain with smooth boundary. Suppose 
that |uo| < 1 and |g| < 1 a.e. in ft and ,/ p '^(uo) < oo. Then the regularized flow 
( 1.20 1 -( 1.22 1 has a unique global weak solution u e S in the sense of Definition 



Moreover, u e ' S satisfies the following energy inequality 



3.1 



j;i{^ 5 { S ))+ i 



\hdt< J^(u ) 



Vs E [0,T]. 



(3.9) 



Proof. Let n p denote any well-known mollifier (cf. [35.), uo iP := r\ p * Uo and 
:= tip * g denote the mollifications of u and g, respectively. Let u e p ' 5 denote the 



classical solution to ( 1.20 1-( 1.22 1 corresponding to the smoothed datum functions Uo iP 
and g p . Hence, u p ,<5 satisfies the energy law (3.8 1 with Uo iP and g p in the place of Uo 
and g. 

we know that |u p ,<5 | satisfies the maximum principle, thus, 



From Lemma 



3.5 



max \u e ' s (x,t)\ < 1. 



Next, since 



lim 

p 



J p°\( u o,p) = J p,aK) < 00 



the energy law for u E p ' S immediately implies that u E p ' S satisfies estimates (3.3)-( [3~6| , 
uniformly in p and 5. 

The remainder of the proof is to extract a convergent subsequence of {u p ,l5 } p >o 



and to pass to the limit as p — > in the weak formulation of (1.201. This can be 



done easily in all terms of (1.201 except the nonlinear term in A p (see ( 1.23 >) . To 
overcome the difficulty, we appeal to Minty's trick or the "decisive monotonicity trick" 
as described in Lemma [2. 1| Since this part of proof is same as that of Theorem 1.5 
of [T5], so we omit it and refer to [T5] for the details. 
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The uniqueness of weak solutions follows from a standard perturbation argument 



and using the fact that is a monotone operator. Finally, the inequality (3.9 1 follows 
from setting p — > in the energy law (3.8 1 for u £ p ' S , and using the lower semicontinuity 
of Jp S x and the L 2 -norm with respect to L? weak convergence. □ 
We conclude this section with some remarks. 

Remark 3.1. (a). Since W ' (Cl) is not a reflexive Banach space, Minty's trick 
would not apply in the case p — 1 without the help of the b p (e)A term in the operator 
At 



- p (see (l.23\). On the other hand, in presence of this term, Minty's trick does 
apply since we now deal with the Sobolev space H 1 ^) instead of W 1,1 ^). In fact, 
if b p (e) = in (1.20), BV solutions are what one can only expect in general for the 



regularized flow ( 1.20 )-( 1.22 ) in the case p = 1 (cf. \2J$ ). 

(b). We also point out that using nonzero parameter b p (e) is not necessary in the 
case 1 <p < oo. For example, the conclusion of Theorem \3. 7| still holds if we replace 
the above b p (e) and p* by 



b p (e) 



if 1 < p < oo, 

e « ifp=l (or>0), 



and 



p 



p if 1 < p < oo , 
2 ifp=\. 



On the other hand, the conclusion of Theorem \3.6\ may not be true any more after 
this modification. 

(c). Theorem 3.1 also holds when CI is a bounded Lipschitz domain, one way to 



prove this assertion is to use the Galerkin method as done in flttf . or to use the finite 
element method to be introduced in Section^ 

4. Passing to the limit as 6 — > 0. The goal of this section is to derive the 
limiting flow of Oob-O^ as S -► for each fixed e > 0. As in [TT1II1II1EI1I1I1I1S], 



the key ideas for passing to the limit are to use the compactness result of Lemma [2^2] 
and to exploit the symmetries of the unit sphere S 1 ™ -1 . Our main result of this section 
is the following existence theorem. 

Theorem 4.1. Letp* := max{2,p} and CI C R m be a bounded Lipschitz domain. 
For 1 < p < oo, suppose that Uo € W l ' p (fi,R n ), |uo| = 1 and |g| < 1 in CI, then 
there exists a map u E G L°°((0,T); W 1 '^ (Cl, R' 1 )) n ^((0, T); L 2 (Cl, R™)) such that 



1 



a.e. in Cli 



(4.1) 



and u e is a weak solution (in the distributional sense, see ( 4.29[ ) ) to the following 
problem 



uf - A e p u £ + A(u £ - g) = ^ A u e inClr, 
B £ p n = on dCl T , 

u £ = u on Cl x {t = 0}, 



(4.2) 
(4.3) 
(4.4) 



whe 



~ b P (s)\Vu £ \ 2 + |Vu £ |r 2 |Vu £ | 2 + A(l - u £ ■ g). 

Moreover, u 6 satisfies the energy inequality 



(4.5) 



■^,a(u £ (s)) + / II u?(t) \\% dt < J| iA (u ) < 4, x (u ) for a.e. s G [0,T], (4.6) 
Jo 
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and the additional estimates 

II ivu e ir 2 vu' 



< C for p — 



P 



div^|| L 2 (i i) < C forl<p<oo, 
div(^Au £ )|| i2(i2) < C forl<p<^, 



1 < p < oo, 



'^(eJllVu'lU-^) < C 
/or some positive e-independent constant C . Here Bp is the n x m matrix 

B £ p := b p (e)Vu £ + |Vu £ |r 2 Vu £ , 

^a(v) :=£?(v) + ^ / W-g\ 2 dx. 



(4.7) 

(4.8) 
(4.9) 

(4.10) 
(4.11) 
(4.12) 



Proof. In light of Remark 3.1 (c) and the density argument, without loss of the 
generality we assume that CI is a bounded smooth domain, Uo and g are smooth 
functions. On noting that |uo| = 1 implies that L 5 (uq) = 0, then, E £ ' S (uq) = E^Uq) 
in (1.161. Since Ep(uo) < E^Uq), Hence, the assumptions on Uo and g ensure (3.2 1 
holds. Since the proof is long, we divide it into three steps. 

Step 1: Extracting a convergent subsequence. Let u £ ' S be t he w eak solution solu- 
tion to ( |1.20 1-( 1.22 I whose existence is established in Theorem 3.7 Since E £ ' S (uq) < 
Ep(u ), then J £ ' x (u ) is uniformly bounded with respect e and 6. Hence, (3.9) im- 
plies that u e - s satisfies the uniform (in both e and S) estimates ( 3.3 1-( 3.6 ) , and the 
maximum principle |u e,l5 | < 1 on CIt- 

By the weak compactness of W 1,p (fl) and Sobolev embedding (cf. [IJ|45]), there 
exists a subsequence of {u e ' <5 } ( 5 > o (still denoted by the same notation) and a map 
u £ G L oo ((0,:r);W 1 ' p *(n,R' 1 )) n ,ff 1 ((0,:r);Z 2 (fi,R n )) such that as o~-> 



,M 



,e,c5 



.Ml 



1 



weakly* in L°°((0, T); W 1,p " (Q, TV 1 )), 
strongly in L 2 ((0, T); L 2 (fL, R")), 
a.e. in fix, 

weakly in L 2 ((0, T);L 2 (il, TV 1 )), 
strongly in L 2 ((0, T); L 2 (f2)). 



It follows immediately from (4.15) and (4.171 that 

|u £ | = 1 a.e. in fiy- 

Step 2: Wedge product technique and passing to the limit. Since |u £ 



(4.13) 
(4.14) 
(4.15) 
(4.16) 
(4.17) 

(4.18) 
< 1 in 



(4.151, that 



fir, an application of the Lebesgue dominated convergence theorem yields, on noting 
strongly in L r ((0,T); L r (n,R n )), VrG[l,oo). (4.19) 



Let 



ce.S 



^(l-|u e 'T)u s ' -A(u' 



,e,<5 
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it follows from the estimate \u e ' 5 \ < 1 and the inequality |u E >*| < |(l + \u £ ' S \ 2 ) that 



o Jn 



fl - lu^l 2 ) • \u E ' s \dxdt< 



o Jn 



fl - |u £ ^| 2 ) 2 dxdt 



ijf ^(i-|u £ 'T)> £ 'T^. 



(3.4) immediately implies that the first term on the right-hand side is uniformly 
bounded in S. Testing equation ( 1 .20 1 by u E ' S and using estimates (|3.3| , ( 3.5 ) and ( 3.6 1 



we conclude that the second term on the right-hand side is also uniformly bounded 
in 5. Hence, i £ ' S is uniformly bounded with respect to 5 in L 1 ((0,T); L 1 (ri,R™)). By 
Lemma [2~2l we have that 



Vu E ^-^Vu e strongly in L 9 ((0, T); L«(ft, R nxm )), Vq€[l,p*). (4.20) 



This and (3.5 1 imply that 



(4.21) 



weakly* in L°°((0, T); L p '(f2, R" xm )) with p' = ^ if p ^ 1 and weakly* in L°°((0, T); 
i°°(fi,R" xm )) if p= 1. 

Next, taking the wedge product of (1.20) with u s ' S yields 



uf A u e ' s - div(B e / A u e ' 4 ) — Ag A u e ' S = 0, 



(4.22) 



where Bt' S is defined in (1.23). 



Testing g22| with any w E L°°((0, T); (f2,R n )) we get 

J {(ufAu^,w) + (8fAu^Vw)-A(gAu E ' S ,w)}^ = 0. (4.23) 



It follows from setting S -> in ( |4.23[ ) and using ( |4.16[ ), ( |4.19| and ( |4.2l[ ) that 
^ |(uf Au e ,w) + (B e p Au e , Vw) - A (g Au e ,w)|rfi = 0, 



(4.24) 



where Bp is given by (4.11 1. 
Note that (14.181) implies 



(4.25) 



uf u e = 0, (S£) J u e = a.e. in Q T - 
This in turn yields the following identity 

J {(uf,uV) + (^,V(u^))+A(u £ -g,u^)}di = ^ (l4, x ,<p)dt (4.26) 

for any <p G i 00 ^) D Z°°((Q,T); W^'P^fi)), where /i 6 ^ is defined by ( gj |. 
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Finally, for any v <E [C 1 (flT)) n , taking w = u E A v in ( |4.24 |, ip = u £ ■ v in ( |4.26| |, 
and using the formula a ■ (b A c) = (a A b) • c yield 

J {(uf , u e A (u £ A v)) + (B e p , V(u £ A (u e A v))) - A(g, u £ A (u E A v)) } dt = 0, 

(4.27) 

(4.28) 



J {(u?,u e (u^v)) + (^,V(u £ (u £ .v))) 

+ A(u e - g, u e (u £ • v)) } dt = £(& x vl b ,v) 



</f. 



Subtracting (4.27) from (4.28), and using the identity 



v = (u e • v) u £ — u e A (u £ A v) 



we obtain that 

i-T 



J {(ul,v) + (B s p ,Vv)+X{u s -g,v)]dt = J (,^ A u £ ,v)dt 



(4.29) 



for any v £ [C^^t)]™- This is equivalent to saying that u e is a weak solution (in the 
distributional sense) to (4.2 1-( 4.4 1 . 



Step 3: Wrapping up. We conclude the proof by showing the estimates (4.6 1 



(4.101. First, (4.6 1 follows immediately from letting S — > in (3.9 1, appealing to 
Fatou's lemma and the lower semicontinuity of L 2 - and L p -norm with respect to L 2 - 
and L p -weak convergence. We emphasize again that this is possible in the case p = 1, 
because the uniform (in S) estimate (3.6 1 implies that u e S L°°((0, T); 7? 1 (S7, R")). 
(4.7 1 and (4.101 are direct consequences of (4.6 1. Finally, the bounds (14. 81) and (4.9 1 



follow immediately from (4.291 and (4.241, respectively. Hence the proof is complete 
□ 

Remark 4.1 

in general when p 



Ifbp(s) = is used in the regularization, the solutions to (4.1 1 -(4.4 1 
are only expected to belong to L°°((Q,T); [BV(f2)] 



1. 



5. Passing to the limit as s —> 0: the case 1 < p < oo. In this section, we 
shall pass to the limit as e — ► in (4.1 1-(4.4 ) and show that the limit map is a weak 



solution to (1.8 1- (1411. Since the analysis and techniques for passing the limit for 



1 < p < oo and p = 1 are quite different, we shall first consider the case 1 < p < oo in 
this section and leave the case p — 1 to the next section. We begin with a definition 
of weak solutions to ( |1.8| -( [l~Tl| in the case 1 < p < oo. 

Definition 5.1. For 1 < p < oo, a map u : Ht — > R" is called a global weak 
solution to (fT8j)-( jl41| ) if 

(i) u € L°°((0, T); W x ' p (n, R' 1 )) n H l ((Q, T);L 2 {tt, R")) ; 

(ii) |u| = 1 a.e. on fix, 

(iii) u satisfies ( |1.8[ )-( i~TT| ) m f/ie distributional sense. 

Our main result of this section is the following existence theorem. 
Theorem 5.2. Let 1 < p < oo, suppose that the assumptions on Uq and g in 
Theorem 4-1 still hold. Then problem ( |1.8[ )-( [l~TT| ) has a weak solution u in the sense 



of Definition\5. l\ Moreover, u satisfies the energy inequality 



J p ,\(u(s)) + f || u t (i) \\% dt < J P: a(u ) for a.e. s e [0,T], 



(5.1) 
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where J Pt \ is defined by (1.7 1. 

Proof. We divide the proof into three steps. 



Step 1: Extracting a convergent subsequence. Let u £ denote the solution of (4.1 1- 



(4.4) constructed in Theorem |4.1| From ( |4.1[ ), (4.6 1— (4.10 1, the weak compactness of 
W l p {Sl) and Sobolev embedding (cf. [Tl 145] ) . there exists a subsequence of {u £ } e>0 
(still denoted by the same notation) and a map u <E L°°((0,T);W 1 < p (il,B. n )) H 
^((O.TjjL^O.R")) such that as e -> 



u £ — u weakly* in L°°((0, T); W 1,p (n, R n )), 
strongly in L 2 ((0, T);L P (D,, R")), 
a.e. in r2y, 



& P (e)Vu e 







weakly in L 2 ((0, T); L 2 (ft, R" xm )), 



u£ — * ut weakly in L 2 ((0, T); L 2 (ft, R")). 



It follows immediately from ( 5.4 1 and (4.181 that 

|u| = 1 a.e. in fl?- 



(5.2) 
(5.3) 
(5.4) 
(5.5) 
(5.6) 

(5.7) 



Step 2: Passing to the limit. First, by (5.4 1 and the Lebesgue dominated conver- 
gence theorem we have that 



u ££ -^u strongly in L r ((0, T); 27(ft, R")), VrG[l,oo). 



(5.8) 



Next, let f e := /i £ A u £ - A(u £ - g). Clearly, f £ € L 1 {(0,T); L 1 ^; R")) and is 
uniformly bounded, on noting (4.5 1-(4.6 1 . By Lemma 2.3 we get 



Vu £ ^Vu strongly in L 9 ((0, T); ^(ft, R" xm )), Vge[l,p), (5.9) 



which, (4.7 1 and (5.51 imply that 

|Vu £ |P~ 2 Vu £ ^ |Vu| p - 2 Vu weakly* in L°°((0,T);L p '(fi,R" xm )), (5.10) 



B £ p ^5 B p := |Vu| p ~ 2 Vu weakly in L 2 ((0, T); L p * (fi, R" xm )), 
where p' = and p'„ := min{2,p'}. 



(5.11) 



It then follows from taking e — > in (4.241 and using (5.6|, (5.8), and (5.11 1 that 

J |(ui Au,w) + (S p Au,Vw) - A(g Au,w)Jdi = (5.12) 

for any w G C 1 ^). Since B p € L°°((0, T);L p '(fl, R" xm )), by the standard density 
argument one can show that fl5~12| also holds for all w G L°°((0, T); W 1,P (Q, R")) n 
L°°(nr). 

Since |u| = 1 a.e. in f2r, there holds the following identity, which is analogous to 



(4.26), 



J |(u t ,u^) + (B p ,V(u<p)) + A(u-g,u</5)|dt = J (n pA ,<f)dt (5.13) 



for any tp G L°°(f2r) D Z°°((0, T); W 1,p (fi,R n )); where ^ P , A is defined by p§. 
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Finally, for any v € [C^f^)]" on choosing w = u A v in (5.121 and cp = u • v 



(5.131, subtracting the resulting equations and using the identity 



v = (u • v) u — u A (u A v); 



we obtain that 

i-T 



{(u t ,v) + (Bp.Vv) + A(u-g,v)}dt = ^ (/i P)A u,v)di (5.14) 



for any v g [C (f2r)]"- Hence, u is a weak solution to ( 1.8 1-( 1.11 1 



Step 3: Wrapping up. We conclude the proof by showing the energy inequality 



(5.1l. First, notice that (4.6 1 implies that 



{Vu £ ( S )P 



a ^P 



-2 



uf (*) < 4,x (uo) VsG[0,T]. 



Then (5.1) follows from taking e 



(5.15) 

in (5.151, using Fatou's lemma and the lower 
scmicontinuity of the L p -norm with respect to L p -weak convergence. The proof is 
complete. □ 

Remark 5.1. (a). We remark that it was proved in \15\ \32[ \4 l^j that weak 
solutions to ( |1.8[ )-( [lTil~| ) are not unique in general. 

(b) . Although the above proof is carried out for any a > in the definition 
of b p (e) (cf. ( 1.15[ > ), a should be chosen large enough so that the error due to the 
perturbation term b p (s)A is much smaller than the error due to other regularization 
terms in numerical simulations. 

(c) . The existence result of Theorem |5.£| is established under the assumption 
u € W /1,p *(f2,R") with p* = max{p, 2}. This condition can be weakened to u € 
W 1,P {Q,, R") in the case 1 < p < 2 by a smoothing technique. 

6. Passing to the limit as s — > 0: the case p = 1. In this section, we consider 
the case p = 1 and establish the existence of global weak solutions for the 1-harmonic 
map heat flow ( |1.8| -( fl~TT| ) by passing to the limit as e — > in (4.2 )-(4.4 1. There are 



two main difficulties which prevent one to repeat the analysis and techniques of the 
previous section. Firstly, the compactness result of Lemma |2.3| does not hold any 
more when p = 1. Secondly, since the sequence {u 6 } e> o is uniformly bounded only 
in L°°((0,T); [W 1 ' 1 ^) n L°°(Cl)] n ), and W^'^R") is not a reflexive Banach space, 
hence, the limiting map u now belongs to L°°((0, T); [BV(to) H L°°(Q,)] n ), i.e., u(t) is 
only a map of bounded variation. As expected, these two difficulties make the passage 
to the limit as e — > become more difficult and delicate. 

6.1. Technical tools and lemmas. In this subsection, we shall cite some tech- 
nical tools and lemmas in order to develop a weak solution concept to be given in the 
next subsection for the 1-harmonic map heat flow. Specially, we need the pairings 
A ■ Dv and A A Dw between a tensor field A and a BV- vector field v, which was 
developed in [53] as a generalization of the pairing b • Dv between a vector field b and 
a JSIZ-function v developed in [4l 110). 

We recall from |23j the space of divergence-L q tensors 



y(n) q := {Ae L°°(n;R nxm ); divA & L 9 (ft;R™)} for 1 < q < oo, (6.1) 



and that A ■ Dv and A A Dv are defined as 
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Definition 6.1. For any A £ y(fl)i and v e [BV(tt) n L°° (fl)] n , we define 
A ■ Dv and A A Dv to be the Junctionals on C^°(f2) and [Cq 3 {Q.)] n , respectively by 

(A-Dv,ip) := - [ (A T v) ■ Vtpdx- [ (divA-v)i/}dx V^eC^(O), (6.2) 
Jn Jn 

(AADv,Vf):=- [ (AAv) -Vwdx- / (div.4 A v) • w dx Vw G [C °° (fi)] n ; (6.3) 
Jn Jn 

where A T stands for the matrix transpose of A and the notation ( 2.4 \ is used in ( 6.3 1 . 

We now list some properties of the pairings A ■ -Dv and A A Dv, and refer to 
Section 2 of [53] for their proofs. The first lemma declares that A ■ Dv and A A Dv 
are Radon measures in Q. 

Lemma 6.2. For any Borel set E C Q, there hold 

\(A-Dv,$)\<max\il>(x)\.\\A\\ La o iEiB » xm) ■ \Dv\(E) W> G C (E), 

| (A A Dv, w)|< max |w(x)| • || A || L ~(£,R- m) ■ \Dv\(E) Vw G [C (E)] n . 

Hence, it follows from the Riesz Theorem (cf. Theorem 1.54 of Jt$) that both func- 
tionals A ■ Dv and A A Dv are Radon measures in f2. 

Corollary 6.3. The measures A ■ -Dv, \A ■ Dv|, A A -Dv, and \A A Dv| all are 
absolutely continuous with respect to the measure |-Dv| in fl. Moreover, there hold 
inequalities 

\(A-Dv)(E)\ < \A-Dv\(E) < \\A\\ L ~ {E ,, R ^ m) ■ \Dv\(E), (6.4) 
\(AADv)(E)\ < \AADv\(E) < \\A lU-^.R-Km) • \Dv\(E) (6.5) 

for all Borel sets E and for all open sets E' such that E c E' c 0. 

Hence, by the Radon-Nikodym Theorem (cf. Theorem 1.28 of J3jj), there exist 
\D~v\-measurable functions O := 0(-4, Dv,x) : f2 — * R, and A := A(A,Dv,x) : f2 — * 
R" such that 

(A-Dv)(E) = [ Qd\Dv\, and || 6 |U=c (n ,|flv|) < II A |U~ ( n,R«* m) ; (6-6) 

JE 

(AADv)(E)= [ Ad\Dv\, and || A ||^ ( o,p v |) < \\A \\ L <~( Q ,nn*™) (6.7) 

JE 

for all Borel sets E d CI. 

The second lemma declares that every A G 3^(^)i has a well-behaved 
traction An on the boundary of a Lipschitz domain Q. 

Lemma 6.4. Let CI be a bounded domain with a Lipschitz continuous boundary 
dfl in R m , then there exists a linear operator (3 : y(Cl)i — > L°°(dQ; R") such that 

||/3(-4)IU-(8n,K») < MIU-(f2,R"x m) , (6.8) 



(A v)an = / v(x) dW" 1 - 1 Vv G [BV(fl) n L°°(Q)] n , (6.9) 

i9(-4)(x)=-4(aj)n(aj) Vie G <9f2, .A G C 1 (?T, R nxm ). (6.10) 



Remark 6.1. Since (3(A) is a weakly defined traction of A on dfl, hence we shall 
use An to denote (3(A) in the rest of this section. 
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The third lemma declares that the following hold. 

Lemma 6.5. Let £1 be a bounded domain with a Lipschitz continuous boundary 
dtt in R m , then for any A G y(Sl)i and v G [BV{Q) n L°°(Cl)] n there hold identities 

[ dWA-vdx+(A-Dv)(n) = [ An-vdH™' 1 , (6.11) 

I dWAAvdx+(AADv)(fl) = I AnfwdH" 1 - 1 . (6.12) 
Jn Jan 

The fourth and fifth lemmas states continuity results for the measure A ■ Dv and 
A A Dv with respect to A and v, respectively. 

Lemma 6.6. Let Aj,A G 3^(^)i and suppose that 

Aj — ► A weakly* in L°°(E), 
div*4j — ► div.4 weakly in L 1 (E) 

for all open sets E CC ft. Then for all v G [BV(Q) H L°°(Q)] n the following hold 

Aj ■ Dv — ► A ■ Dv weakly* in M(ft), (6.13) 

Aj A Dv — ► A A Dv weakly* in [M(Q)] n , (6.14) 

and 

e(Aj,Dv,-) — >Q(A,Dv,-) weakly* in L°°(E)for all E CC fl, (6.15) 

A(A v Dv,-) — >A(A,Dv,-) weakly* in[L°°{E)] n for all E CCO. (6.16) 

Here "E CC Q" means that E is compactly contained in f2; that is, E C E C £1 and 
E is compact. 

Lemma 6.7. Let A G y(to)i and v G [BV(tt) H i°°(fi)] n . Suppose that {vj} C 
[C°°(fi) n flV(fi)]" s*ric% converges to v (cf. Definition 3.14 of J^). Then 

A ■ Dvj — > A ■ Dv weakly* in M(ft), (6.17) 
A A Dvj — ► A A Dv weakly* in [M(Q)] n . (6.18) 

Moreover, 



A-Dvjdx — > / A-Dv, (6.19) 
Jn 

A A Dvj dx — > I A A Dv. (6.20) 



Jn 

The sixth lemma gives the precise representations for the density functions O and 
A defined in Corollary |6.3| 

Lemma 6.8. (i) Lf A G y(n) 1 n C*(fi,R nxm ) ond v G [BV{Cl) n L°°(fi)]™, t/ien 
i/iere ZioW 

e(A,Dv,x)=A(x)-j^-(x), A(A,Dv,x)=A(x)A—{x) (6.21) 
|Dvl — a.e. in 
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(ii) IfAe y(fl)i and v e [BV(to) n L°°(to)] n , then there hold 

Q(A,Dv,x)=A(x)-^ 1 (x), A(A,Dy,x) = A{x)A^r(x) (6.22) 
\Dv\ \Dv\ 

|Dv| a — a.e. in to. Where denotes the density function of the measure Dv with 
respect to the measure \Dv\, and |Z?v| a denotes the absolute continuous part of the 
measure \Dv\ with respect to the Lebesgue measure C n . 

Next, we recall from |23j the space of divergence-measure tensors 

VT(Q) := {At L°°(ft,R nxm ); div.4 e [M(fl)] n }, (6.23) 
and briefly discuss two of its important properties. Firstly, like in the case of the 



space of the divergence-L 1 tensors y(to)i, Definition 6.1 is still valid for A € T>T(to) 
andv G [BV(to)r\L°°(to)nC(to)] n (cf. g]). Secondly, there is a well-behaved traction 
An for every A € VT(to). 

Lemma 6.9. Let £1 be a bounded domain with Lipschitz continuous boundary dto 
in R m , then there exists a linear operator a : T>T(to) — > L°° (9f2;R™) such that 

II a(A) |U» (a n,R») < II A |U»(fi,Rn*m), (6-24) 

{A,v)ga= I a(A)(x)v{x)dH m - 1 Vv e [BV(to) n L°°(fi) n C(to)] n , (6.25) 
Jan 

a{A){x) =A(x)n(x) Vx G dto, A e C\to, R" xm ). (6.26) 

Moreover, for any A G VT(to) and v G [BV{to) n L°°{Q) n C(to)] n , let An := a(A) 
on dto, then there hold the following Green's formulas 

(divA-v)(to) + (A-Dv)(to) = [ An-vdH™- 1 , (6.27) 

{divAAv){to) + {AADv){to) = f AnAvdH™- 1 . (6.28) 

Jan 

Thirdly, there holds the following product rule . 

Lemma 6.10. For any A G T>T(Q) and v G [BV(to) H L°°(Q)] n , the identities 



div(^v) = (div.4) • v + A- Dv, (6.29) 
div(„4 Av) = (div„4) AvliADv. (6.30) 

hold in the sense of Radon measures in to. Where v denotes the limit of a mollified 
sequence for v through a positive symmetric mollifier, A ■ Dv (resp. AADv) is 
a Radon measure which is absolutely continuous with respect to the measure |-Dv|, 
and whose absolutely continuous part (A-Dv) a (resp. (A/\Dv) a ) with respect to 
the Lebesgue measure C m in to coincides with A - (Vv)" (resp. A A (Vv) a j almost 
everywhere in to, that is, (A ■ Dw) a = A ■ (Vv) a (resp. (A A Dv) a = A A (Vv)"j, 
C m — a.e. in to. 

Remark 6.2. ( |6.29[ ) and (6.301 hold without all the overbars if either v G 

[BV(to) n L°°(ft) n C(to)} n orAe y(fi)i. 

Remark 6.3. Lt should be note that the results of Lemmas \6.6§6.8\ also hold for 
the tensor fields in VT{to) and the vector fields in [BV(to) n L°°(to) n C(to)] n . 
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6.2. Existence of weak solutions of 1-harmonic map heat flow. Through- 
out the rest of this paper we let Bi (R nxm ) denote the unit ball in the Euclidean space 

Rnxm. that jg^ 

m n i 

Bi(R „x m) = j^gRnxm. 1^ ;= < l}. 

In addition, let 

Vu £ 

A : = Ml,A u ) O 



Vu e 



Hence B\ = 0l (e)Vu £ + B £ (cf. ( jiH] ). 

We now give a definition of weak solutions to ( |1.8p - p~TT| ) in the case p = 1. 
Definition 6.11. Forp = 1, a map u : f2-r — ► R™ is called a global weak solution 



to ( |l.8[ )-( |l.ll[ ) if there exists a tensor (or matrix-valued function) B such that 

(i) u g i°°((0, T); n L°°(fi)] n ) n ^((0, T);L 2 (n, R")), 

(ii) |u| = 1 £ m+1 -a.e. zn fi r , 

(hi) 6€ L°°((0, T); L°°(!], B^R"™))) n L 2 ((0, T); PT(fi)), 

(iv) u and ,8 satisfy B hue L 2 ((0, T); y(f2) 2 ), B T u = Q, and B ■ (Du) a = |flu| a 
£ m+1 -a.e. infir. 

(v) fin = on 9J7t in the sense of Lemma \6~. 

(vi) £/iere ZioZds i/ie identity 



J J |u t -v + B- Vv + A(u-g) ■v^dxdt= J (^Jvda x )<ll 



for any v £ C^fir)- Where o\ denotes the vector-valued Radon measure 

cr A = (fiAu)Aflu+A(l-g-u)u. (6.31) 

Moreover, the Radon measure (B A u) A flu is absolutely continuous with 
respect to the measure \Du\, and for C 1 — a.e. t £ (0,T) there exists a 
function $(B, flu, x,i) : CI — > R such that 

((BAu)Aflu)(£)= / $(fi A u, flu, x,i)d|flu| for all Borel sets E C n, 

J E 

ll*(*)IU»(n,|23u|)<||fl(*)IU»( n , R «««) /or/: 1 a.e. t€(0,T), 

$(B Au,Du,x,t) = (S(t) Au) A — ^ |flu| a a.e. m 0, £ x a.e. t £ (0,T), 

((B A u) A flu)" = (B A u) A (flu) a = |Au| a u C m+1 a.e. in fi x . 

Remark 6.4. (a). The tensor field B extends Tjj^i as a calibration across the 
gulfs. 

(b). ifu(t) g w{£(p) and Bit) £ y{n) x , using (ii) and (iv), and the identity 

(a A b) A c = (a • c)b — (b • c)a we get that (B A u) A flu = |Vu| u. Hence, (6.31 ) can 
be rewritten as 

<7 X = |Vu| u+ A(l - g • u)u. (6.32) 

Thus, (6.31) is a weak form of ( |6.32 1 since |flu| u may not be defined for u(t) g 
[BV(n) n L°°(n)} n and B{t) £ VTjtf). 
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Our main result of this section is the following existence theorem. 
Theorem 6.12. Let p = I, suppose that u G H x (fi, R^) a nd g G L 2 (fi,R n ) 
wii/i | Up | = 1 and |g| < 1 a.e. in Q. Then, problem ( |1 .8[ > -( [TTTT] ) /ias a global weak 



solution u in the sense of Definition 6.11 Moreover, u satisfies the energy inequality 

(6.33) 



whe 



h(u(s)) + f || u t (t) ||| 2 d< < J A (u ) /or a.e. s G [0,T], 
Jo 

/ A (u) := \Du\(Sl) + ^J |u-g| 2 dx. 



(6.34) 



Proof. The proof is divided into four steps. 

Step 1: Extracting a convergent subsequence and passing to the limit. Let u E 



denote the solution of (4.1 1-( |4.4| constructed in Theorem 4.1 It follows from (4.1l, 
{46]), ([17]), and ([18} that {u e } £>0 is uniformly bounded in L°°((0,T); [W 1A (n) n 
L°°(n)} n ) n J ff 1 ((0,? 1 );L 2 (fi,R™)), {£ £ } e > in L°°((0, T);L°°(Cl, Bi(R" xm ))) and 
{div£f} e>0 in L 2 ((0, T); L 1 (0; R n )), and {o-{} £ >0 is uniformly bounded in L°°((0,T); 
i 1 (n,R n )). Since L 1 ^) C 7W(fi) and FF 1 ' 1 ^) C W(ft), by the weak com- 
pactness of A4(£l) and -BF(f2) (cf. [5]) we have that there exist subsequences of 
{u e } e> o, {B £ } e> o, {2?f} £ >o, and {<J £ x } E> o (still denoted by the same notation), re- 
spectively, and maps u G L°°((0,T); [BV(Sl) n L°°(f2)] n ) ni/ 1 ((0, T); L 2 (f2, R")), 
6 G i oo ((0,r); J L oo (O, J Bi(R" xm ))), i/ G L 2 ((0, T); [X(0)j"), and <r A G L°°((0,T); 
[yM(f2)] n ), respectively, such that as e — > 





— > u 


weakly* in i°° ((0, T) ; n (fi)] ") , 


(6.35) 






strongly in L 2 ((0, T);£ x (f2, R")) 


(6.36) 






a.e. in fix, 


(6.37) 


6i(e)Vu e - 


-> 


weakly in L 2 ((0, T); L 2 (fi, R™ xm )), 


(6.38) 


uf- 


-> u t 


weakly in L 2 ((0, T); L 2 (fi, R™)), 


(6.39) 


divSf - 


— * 


weakly* in i 2 ((0,T); [M(Q)] n ), 


(6.40) 




-» 6 


weakly* in i°°((0, T); L°°(r!, B 1 (R" xm ))), 


(6.41) 






weakly* in L°°((0,T); [M(O)]™). 


(6.42) 



It follows immediately from (6.371 and ( 4.18[ ) that 

|u| = 1 a.e. in fi<r, 



(6.43) 



and an application of the Lebesgue dominated convergence theorem yields that 

u £ ^5u strongly in L r ((0, T); L r (ft, R")), Wg[1,oo). (6.44) 

Now, taking e 
[C 1 ^)]™ and w G i oo ((0,l7[ff I (n;R n )7 

l-T 



in ( |4.24| ), and ( |4.29| (with p = I) we have for any v G 
0,T)~e%};R™)) 

y j(u t A u, w) + (S A u, Vw) - A (g A u, w)| 



dt 



{( 



B,Vv) + A(u- 



I.v)} 



= 0, 
vdcr A 



(ft. 



(6.45) 
(6.46) 
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Step 2: Identifying v and g\. Firstly, it follows from the identity (B £ ) T u £ = 
(cf. flOS) ), ( |6lT| ), and ( pHIf that 

6 T u = a.e. in ft T . (6.47) 

Secondly, for any v <E [C 1 (Q T )] n , it follows from ( |6^40| ), and ( [6~4"l| ) that 

/ v ■ dv = lim / v • divSf efedt = — lim / Vv • B\ dxdt = — I Vv • Bdxdt. 

Hence, divB exists and 

div£=^, therefore B e L 2 ((0, T); T>T(Cl)). (6.48) 



It then follows from (6.11 1 that 



Bn = 



Thirdly, notice that (6.451 immediately implies that 

div(SAu) S L 2 (n T ,TL n ), hence B A u E L 2 ((0, T); y(tl) 2 ). 



(6.49) 



Let {u p } denote the smooth approximation sequence of u as constructred in 
Theorem 3.9 of [3]. For any v € [C^(Ct T )] n , set w = u p A v in ( |6.45[ ) we get 

J |(u t ,u A (up A v)) + (B Au, V(up A v)) - A(g,u A (u p A v))jcft = 0. (6.50) 

It follows from (6.43), the convergence property of u p (cf. 3 ), and the identity 
u A (u A v) = u(u • v) — v that 

lim(u t ,uA (u p Av)) = (u t ,uA (uAv)) = (u t ,u(u- v) - v)) = -(u t ,v), 

lim (g, u A (u p A v)) = (g, u A (u A v)) = (u - g, v) - ((1 - g • u)u, v) . 



It follows from Theorem 3.9 of [3], Lemma 6.7 and the identity 

(B A u, V(u p A v)) = ((B A u) A u p , Vv) + ((B A u) A Vup, v) 

that 

lim (B A u, V(u p A v)) = ({B A u) A u, Vv) + ((B A u) A Du, v) 
= -(B,Vv) + ((B Au) A Du, v). 



Here we have used the fact that (6Au)Au = — B in the light of (6.43) and (6.471, 
and the measure (BAu)A Du is defined by (6.3 ) with A = B A u and v = u. 

Finally, substituting the above three equations into (6.501 and multiplying the 
equation by (—1) we get 

J {(u t ,v) + (B,Vv) +A(u-g,v)}df=^ ((SAu)ADu (6.51) 

+A(1 - g • u)u, v) dt 
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for any v € [C^(Q T )] n - This and ( E46) imply that 

cr A = (# A u) A Du + A(l - g • u)u. 

S'iep 5: Identifying B. First, since {B e ■ Vu E } is uniformly bounded in L 2 ((0,T); 
L 1 (fi)), then there exists a subsequence (still denoted by the same notation) and 
H £ L 2 ((0,T);.M(ft)) such that 



B E • Vu E 



fj. weakly* in L 2 ((0,T); M(n)). 



(6.52) 



Let u e p and u p denote mollified sequences for u e and u, respectively, through a 
positive symmetric mollificr. For any open set E C f2 we have 



lim 



B 6 ■ Vu e dxdt = lim lim 

, T 

lim lim 

p^O e^O 
T 







B e ■ Du p dxdt 



B 6 ■ Du E p dxdt 



(6.53) 



= lim / (B-Du p )(E)dt (by (IfUl), d6l44l) 



= / B~Du(E)dt, (by (Egg)) 
Jo 

where $ • Du is defined in Lemma 

Hence, it follows from (6.521, (6.53) and Lemma 6.10 that 



fi = B Du« \Du\. 

We refer the reader to Definition 1.24 of [3] for the notation "<< r 
On the other hand, a direct calculation yields that 

|Vu £ i 2 



(6.54) 



B e • Vu e 



> |Vu £ 



Setting e — > and by the lower semicontinuity of the BV-^aoxm (cf. pJj) we obtain 

fi » \Du\, 



which together with (6.54) and Lemma 6.10 yield that 

\Du\ = n = B- Du, 

hence, 



(6.55) 



|-Du| a = (B ■ Du) a = B ■ (Du) a . (6.56) 
Finally, we note that all other properties of B and the measure (B A u) A Du 



listed in (vi) of Definition 
Ol [OplOl and Corollary 



6.11 are immediate consequences of (6.551 and Lemmas 



6.3 



Step ^: Finishing up. We now conclude the proof of the theorem by showing the 



energy inequality (6.331. First, notice that (4.6) implies that for a.e. s 6 [0,T] 



[ s \\um\\h dt<ji x (u ) 

Jo 



(6.57) 



VuoHi* +ai(e)|n|+I A (uo). 
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Then (6.33) follows from taking e — > in (6.571, using the lower semicontinuity of the 
£?V^-seminorm and the L 2 -norm with respect to L 2 -weak convergence. The proof is 
complete. □ 

Remark 6.5. (a). Since weak solutions to ( |1.8| -( fETT) are not unique in general 
for 1 < p < oo (cf. \15[ \3Sl \4^V , we expect that this nonuniqueness also holds for the 
case p = 1. 

is proved under the assumption u G H (fi, R"). 



(b). The existence of Theorem 



6.12 



This assumption can be weaken to Uo G [BV(£l) PI L°°(fl)] n using a smoothing tech- 
nique. 

7. Fully discrete finite element approximations. 

7.1. Formulation of fully discrete finite element methods. For ease of 
exposition, we assume f2 is a polytope in this section. Let T/, be a quasi-uniform 
"triangulation" of the domain of mesh size < h < 1 and = [J KeTh K (K e Th 
are tetrahedrons in the case m — 3). Let J T :— {tk}^ =0 be a uniform partition of 
[0,T] with mesh size r := ^, and 9 t v fe := (v fc — v k ~ 1 )/r. For an integer r > 1, let 
P r (K) denote the space of polynomials of degree less than or equal to r on K. We 
introduce the finite element space 

V' 1 = {w h e C(n,TL n ) n H 1 (f2,R n ); v h \ K e [P r (K)] n , \/K e %}. 



Notice that the density function F defined in (1.171 is not a convex function. On 



the other hand, there exist two convex functions W+ and W- such that 



F(v) = W+(v) - W-(v). 
and W_(v) 



■i 



2 



(7.1) 

i. Clearly, the above 



One such an example is W + (v) 
decomposition is not unique. 

We are now ready to introduce our fully discrete finite clement discretizations for 
the initial boundary value problem (1.20 1-( 1.22 1. Find e V' 1 for k = 1,2,- •• , L 
such that 



(dtulvh) + (BlVv h )+\{u k h -g,v h ) 
B^=[6 p ( £ ) + |Vu^r 2 ]Vu^ 



(7.2) 



WKuJJ-^.vfc) VveV", 



(7.3) 



with some starting value u° S to be specified later. Note that for notational 
brevity we have omitted the indices e, S and p on u k and B\ 

For each k, equation (7.2 1 is a nonlinear equation in u^. Hence, the above nu- 
merical method is an implicit scheme, its well-posedness is ensured by the following 
theorem. 

Theorem 7.1. For each fixed k > 1, suppose that u^ _1 £ ~V h is known, then 
there exists a unique solution u k E V' 1 to (7.2 1 -(7.3). Moreover, {u^}£_ satisfies the 
following energy estimate 



2 ^ 

fc=i 



^K)<W forl<i<L. 



re, 5 



(7.4) 



-ffere ,7^ is defined by (1.191 
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Proof. For each fixed k > 1, it is easy to check that ( 7.2 )-( 7.3 1 is the Euler- 
Lagrange equation of the following functional over V h 



G fc (v) 



2r' 



iJ- l | a + ^|Vv| a + i|Vv|f + ^|v 



(7.5) 



+ ^W+(v) } dx - i y IT (u*- 1 ) -vdx. 



Since is a convex, coercive and differentiable functional, then it has a unique 
minimizcr ujj £ V 71 (cf. 



), hence, (7.2)-(7.3) has a unique solution. 

7T^ 



Since is the minimizer of Gk over V , we have that 

G fc (u£) < Gkiu^ 1 ). 
It follows from the convexity of W— that 

w^.(u*- x )(uj; - uf- 1 ) < W.(u k h ) - W-iut 1 ). 



(7.6) 



This and (7.6 1 imply that 
1 



< 0. 



The bound (7.4 1 then follows from applying the summation operator T^2 k=1 (1 < £ < 
L) to the last inequality. Hence the proof is complete. □ 

7.2. Convergence analysis. The goal of this subsection is to show that the 



numerical solution of ( 7.2 )-( 7.3 1 converges to the unique weak solution of (1.201 



(1.221 as h,r — * 0. There are two approaches to reach this goal. The first approach 



3.7 



assumes the e xistence of the solution of (|l.20|-(|l.22|), which in fact has been proved in 
Theorem 



and then proves that u£ converges to that solution. The other approach 



shows the convergence without assuming the existence of the solution of ( 1.20 1-( 1.22 1 



This can be done by a pply ing the energy method and compactness argument used in 



the proof of Theorem 3.7 to the finite element solution {u^}. In the following, we 



shall go with the latter approach since this will also provide an alternative proof for 



Theorem 3.7 as alluded to in Remark 3.1 (c). 

For the fully discrete finite clement solution {u^}, we define its linear interpolation 
in t as follows 



tk - 1 k _i 



(•) Vtefo-i.ffc], l<k<L. (7.7) 



Clearly, TJ e - s < h ' T is continuous in both x and t. 

The main result of this section is the following convergence theorem. 

Theorem 7.2. For 1 < p < oo. suppose that u £ W 1 ^* (ft, R"), |u | = 1 and 
|g| < 1 infl. For each pair of positive numbers (e,S), let XJ 6 < s > h ' T be defined by (7.7 1. 
Then, there exists u £ > s £ L°°(fi T ) such that 



lim II u e 

fe,T->0 



U 



|£«(0 T ) 







Vg £ [l,oo), 



(7, 



provided that 



lim 



u 



1 h \\W 1 -p*(Q) 
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Moreover, u e,<5 solves ( 1.20 )-( 1.22 1 in the sense of Definition 3.1 



Proof. The proof follows the same lines as that of Theorem 1.5 of [IS], where 
the convergence of a general Galerkin approximation was proved for the case p > 2. 
Since the finite element approximation is a special Galerkin approximation, the proof 
of Theorem 1.5 of [12] can easily be adapted to the finite element approximation 
jje,S,h,T f or p -> 2 thanks to the discrete energy estimate (7.4 1 and the facts that 
dtU * = Uf> h ' r and tYLi W d t<\ 



2 

L- 



= IIU 



:,8,h,T 



\L 2 (L 2 



Since the operator — is uniformly elliptic, as a result, the compactness of 



Lemma [2. 2 1 not only holds for p > 2 but also for 1 < p < 2. In addition, note that 
w ,sm,t ig uniformly ( in h an d r ) bounded in L°°((0, T); H^Sl, R n )) for 1 < p < 2. 
Hence, the proof of Theorem 1.5 of [IS] can be adapted with slight modifications to 
prove (7.8 1 for the case 1 < p < 2. □ 

Remark 7.1. Several practical choices o/u° are possible. For instance, both the 
L 2 -projection of Uq and the Clement finite element interpolation of Uq into *V h (cf. 
H4V are qualified candidates for u°. 

An immediate consequence of Theorems |4.1[ |5~2| |6.12| |7.1| and |7.2| is the following 
convergence theorem. 

,S,h- 



Theorem 7.3. Let 1 < p < oo and V £ ' d ' h - T be defined by (|7.7l 
assumptions of Theorems \T7\ \7J% gTTJ \5J^ and\6J^hold. Th 



the 



ten, there exists a sub- 
sequence of {U £ ' S ' h ' T } (still denoted by the same notation) and a weak solution u of 
|D!])-( |1.11[ ) such that 



lim 



lim 

7i,t->0, 



U 







V?g [l,oo). 



(7.9) 
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